NASA Technical Memorandum 87756 


Buckling Characteristics 
of Hypersonic Aircraft 
Wing Tubular Panels 


William L. Ko 

Ames Research Center 

Dry den Flight Research Facility 

Edwards , California 

John L. Shideler 

Langley Research Center 
FI amp ton, Virginia 

Roger A. Fields 
Ames Research Center 
Dryden Flight Research Facility 
Edwards, California 


NASA 

National Aeronautics 
and Space Administration 

Scientific and Technical 
Information Branch 


1986 


Introduction 

When the speed of an aircraft reaches the hy- 
personic range, aerodynamic heating becomes severe. 
Some of the hypersonic flight vehicle structural con- 
cepts that have been advanced use a thermal pro- 
tection system (TPS) to prevent the structure from 
overheating. For example, the Space Shuttle uses a 
TPS designed to limit the structural temperature to 
350°F (a warm structure). A different concept pro- 
posed for future hypersonic aircraft (ref. 1) was an 
aerodynamically acceptable wavy heat shield made 
of heat resistant metal, such as Rene 41, to limit 
the structural temperature to about 1350°F (a hot 
structure). Compression buckling is a major concern 
for hot structures because of the combined effects of 
aerodynamic load, thermal stress, and reduction in 
material moduli (i.e., modulus of elasticity E and 
shear modulus G ). 

Studies of structures for future hypersonic flight 
vehicles have identified advanced structural concepts 
which show promise of having low structural unit 
mass and high buckling strength (refs. 2 to 13). 
Since curved shell sections exhibit high local buck- 
ling strength, most of the structural panel concepts 
investigated used curved surfaces to achieve high 
buckling strength. Two of the hot structural panel 
concepts investigated were beaded panels and circu- 
lar tubular panels (ref. 1). Results from extensive 
buckling studies of Rene 41 beaded panels and alu- 
minum circular tubular panels are reported in refer- 
ence 9 and references 2 to 6, respectively. All the test 
results for tubular panels were obtained from single- 
panel, room-temperature loading tests under 
laboratory conditions. Furthermore, only limited 
buckling data have been reported on Rene 41 non- 
circular tubular panels, which (based on the analysis 
methods of reference 6) are more efficient than circu- 
lar tubular panels for lightly loaded conditions. 

Thus, to characterize the buckling behavior of the 
tubular panels under combined loads and at elevated 
temperatures, five Rene 41 noncircular tubular pan- 
els (fig. 1) were attached to the wing root region 
of the hypersonic wing test structure (HWTS, de- 
scribed in the Test Equipment section). These five 
panels, which replaced beaded wing panels on the 
HWTS (refs. 9, 12, and 13), were exposed to ex- 
tensive nondestructive buckling tests under different 
combined load conditions (axial compression, bend- 
ing under lateral pressure, and shear) at uniform 
temperatures of 70°F, 550°F, and 1000°F. The use of 
uniform elevated test temperatures caused a reduc- 
tion in material moduli (E and G ) and minimized 
thermal stresses due to temperature gradients. Al- 
though they are difficult to calculate and measure, 


thermal stresses will undoubtedly have to be consid- 
ered in the future if hot structures are to be used. 
The buckling loads were estimated through use of 
the force/stiffness (F/S) method of plotting the test 
data (ref. 11), and the results were compared with 
theoretically predicted buckling interaction curves. 
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radius of circular arc of panel tube 
cross section 

ratios of actual stress to critical stress 
for bending, compression, and shear 

output of rosette strain gage 
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panel tube cross section, 2 aR 
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Five identical Rene 41 tubular panels were de- 
signed and fabricated to replace the root-chord wing 
panels of the HWTS. The design of the panels, de- 
scribed in references 2 and 6, used a random search 
optimization routine to determine values of the cross- 
section variables which constitute a minimum mass 
per unit area subject to specified applied load, geo- 
metric, and failure constraints. The panel design 
loads were N x — 800 lbf/in., N xy = 250 lbf/in., and 
p = 0.75 psi at T = 1350°F. The resulting design, 
which was constrained by a minimum skin thickness 
of 0.016 in., is shown in figure 1. Although the av- 
erage thickness of each chemically milled sheet was 
determined to be 0.0168 in., the design thickness of 
0.016 in. was used in the analyses throughout this 
paper. 

Each tubular panel was made of two formed 
Rene 41 alloy sheets seam welded together to form 
five flat regions (double sheets) and four noncircular 
tubular regions (i.e., flattened tubes). The Rene 41 
was procured in a solution annealed (1975°F) condi- 
tion. Prior to welding, the circular arcs in each sheet 
were incrementally brake formed, and the end clo- 
sures were die formed. Doublers were spot welded to 
both sides at each end of the panel to prevent local 
end failure and to reduce excessive deformations due 
to shear loads. After the final weld assembly, the 
panels were aged for 1 hour at 1650°F followed by 
10 hours at 1400°F. Figure 2 shows a photograph of 
one of the fabricated tubular panels. The panel had 
eight attachment points for z-shaped clips to support 
the heat shields, which are described subsequently. A 
detailed description of the panel fabrication process 
is given in references 3 and 6. 

Analysis 

Local instability is, by design, the critical failure 
mode for the tubular panel shown in figure 1. Al- 
though local and general instability are nearly equal 
under some combined load conditions (ref. 2), it is 
likely that local instability would occur at the same 
time, even if general instability were to occur first. 
Consequently, this paper primarily addresses local 
buckling behavior and the equations governing lo- 
cal buckling. (The general buckling equations used 
in the design and analysis of the tubular panel are 
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those identified by Euler (wide column) for compres- 
sion and by Timoshenko for shear. The equations, 
which include the effects of plasticity and bending 
due to an initial imperfection, are given in refs. 2 
and 6.) 

Local Buckling 

For a tubular panel with tubes of completely cir- 
cular cross sections, the equations for local buckling 
(bead crippling) of circular arc elements of the panel 
in compression, bending, and shear may be written 
in the current notation as follows: 

Compression (eq. (14-3) of ref. 2): 

ft \ 119 

fcc = 0 . 738 % E{-) ( 1 ) 

Bending (eq. (12-33) of ref. 2): 

/n 1 * 15 

fcb = 0-77r/ 3 E (— ) (2) 

Shear (eq. (12-34) of ref. 2): 

fcs = mGK a (Jr) 2 (3) 

where 

/ q2 , \ 0 ' 514 

= (4) 

Buckling equations (1) and (2) are valid for the range 
20 < -y < 120, and equation (3) is valid for 

>50 

To apply equation (1) to the noncircular tubular 
panel, a knockdown factor of 0.86 is recommended. 
(See ref. 6, p. 46.) Equations (2) and (3) are applied 
directly to the flattened tubular panel without mod- 
ification. The buckling equation for compression of 
the noncircular tubular panel is then 

/ t \ 119 

/ cc = (0.86)(0.738)r /3 E^-j 
or 

ft \ lA9 

f cc = 0.635% (5) 

Equations (5), (2), and (3) are used to calculate the 
theoretical buckling strains f cc /rj^E , f c b/WdE, and 


fcs/V2G in compression, bending, and shear, respec- 
tively, for use in the force/ stiffness plots of the non- 
destructive buckling data described in a subsequent 
section. 


Compression-Shear Interaction 

The standard interaction equation for buckling 
failure of a panel under combined loads of axial 
compression and shear is (ref. 2) 

Rc + Rg — 1 

where R c and R s are ratios of the actual compres- 
sion and shear stresses in the panel at failure under 
combined loads to the critical stresses in pure axial 
compression and in pure shear, respectively. This 
equation is used for all buckling failure modes. For 
general instability, the stress ratios are defined as 


Rc = 


Nx 

(^Yr)cr 


Rs 


N xy 
(N X y) cr 


For the local buckling mode for 
tubular panel (bead crippling), the 
defined as 


d _ fc , fb 

Rc ~7Z + T„ 


the noncircular 
stress ratios are 


Rs 


A 

fcs 


The stress ratio for local buckling in compression ac- 
counts for coupling between compression and bend- 
ing. This coupling occurs even when zero lateral 
pressure is applied to the panel because an assumed 
initial imperfection of 0.001L provides a moment arm 
by which compression can always produce a bending 
stress. (See ref. 2.) 


Plasticity 

The plasticity correction factors which appear in 
the local buckling equations are defined as (ref. 2, 
p. 31) 

Vl — ^tan = (6) 

V2 — Vsec — Rsec/R (7) 

m = = (%an»?sec) 1/2 (8) 

Through use of a modified Ramberg-Osgood stress- 

strain approximation (see ref. 2), the tangent and 
secant moduli at the buckling stress are 


L'tan 


/cr 

(/cr/ E') + n[0.002(/ cr / f C y) n - 0.00001] 


(9) 
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and 


Esec — 


fcv 

ifcr/E) + 0.002(/ cr // cy ) n - 0.00001 


( 10 ) 


where 


and 


/cr > /pi — / C y (0.005) — 


^tan — Esec — E if / cr < fp\ 


( 11 ) 


The shape factor n in the Ramberg-Osgood stress- 
strain approximations for Rene 41 (eqs. (9), (10), and 
(11)) is taken as n = 25.0 at 70°F, 22.2 at 550°F, and 
18.5 at 1000°F. Figures 3, 4, and 5 show the plots 
of equivalent elastic stress fcr/Vi (i = 1,2,3) as a 
function of actual stress / cr for the three respective 
temperatures 70°F, 550°F, and 1000°F. Values of the 
modulus of elasticity E and the shear modulus G 
for Rene 41 are shown in figure 6 as a function of 
temperature (ref. 14). With the aid of figures 3 to 6, 
the theoretical buckling strains in compression (e c )cr> 
in bending (c&) C r, and in shear 7 cr can be calculated 
from equations (5), (2), and (3) as 


(^c)cr — 


fee 

m E 


(12) 


( e b)cr 


feb 

m E 


7cr — 


fes 

mG 


(13) 

(14) 


The values of (e c ) cr , (e{,)cr, 7cr, V2, V3, E, G, and / cy 
for the different temperatures are given in table 1. 


Force/Stiffness Method 

The purpose of conducting nondestructive buck- 
ling tests instead of destructive buckling tests was 
to avoid the cost associated with destructive tests of 
a large number of panels. In destructive buckling 
tests, only one buckling data point for one load con- 
dition can be generated from each test panel. How- 
ever, through use of the F/S method to predict the 
buckling strength, a wide range of buckling data can 
be generated from each test panel for different load- 
ing and temperature conditions. The F/S method 
was advanced by Jones and Greene (ref. 11) for the 
prediction of general and local buckling strengths of 
structural components whose buckling behavior is 
complex or nonlinear. 

Since local buckling is, by design, the failure 
mode for the tubular panels, the F/S method used in 
this paper is one developed to predict local buckling 
failure. The method uses a plot of F against F/D , 


where F is the applied load and D is a generalized 
strain variable which accounts for axial compression, 
bending, and shear components. The generalized 
strain variable D is given by 



and the predicted local buckling occurs when 

D = D ct = 1 and ?- = ^=F ct (16) 
U yJqv 

The strains e c , e and 7 are measured with strain 
gages, and the buckling strains (e c ) cr , ( € 6)cr, and 
7 cr are calculated from equations (12), (13), and 
(14), respectively. (See table 1.) The exponent m 
in equation (15) was empirically determined to be 
2 for most types of panels including the completely 
circular tubular panels (ref. 5). For the present F/S 
analysis, m = 2. 

Equations (15) and (16) represent a buckling 
strain interaction surface which is the basis for the 
limit strain lines used in the F/S plots. Figure 7 
shows a graphical illustration of the F/S method, 
which requires extrapolation of the curve fitting the 
test data points. The buckling failure load is deter- 
mined from the intersection of the extrapolated curve 
and the limit strain line. The accuracy of buckling- 
failure-load prediction with the F /S method depends 
on (1) the location of strain gages so that they mea- 
sure strain which is sensitive to the impending buck- 
ling mode shape, (2) the distance of extrapolation 
(that is, how close the final test data point is to the 
limit strain line), (3) the accuracy of the curve fitting, 
and (4) the accuracy with which the limit strain line 
itself is determined (e.g., if the critical strains are 
determined experimentally for a specific configura- 
tion, they may be more accurate than if they are 
determined analytically from general equations). 

The extrapolation of the F/S test data points 
to the intersection with the limit strain line was 
accomplished through least-squares fitting of the test 
data through use of the following equation from 
reference 11: 

F = 1 + CiF + C 2 ^ 2 

D C 3 + C 4 F + C5F 2 + C 6 F s 1 ; 

The buckling value of F (the intersection point Eqt) 
was determined by setting D = 1 in equation (17). In 
the present F/S analysis, D was expressed as follows: 
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Room temperature: 


then (N x ) cr and (N xy ) cr can be estimated from 


n SG,+SG, , |SG, -SOjI , /'2|RSG 2 -RSG,|V 

d = ^uT + + ( — J (18) 

Elevated temperatures: 


D = 


(SGi -SGi)-f (SG 2 -SG 2 ) 

2(cc)cr 

f |(SG 1 -SGx)-(SG 2 -SG 2 )1 

2 ^b) cr 


2|(RSG 2 - rsg 2 ) - (RSG 3 - RSG 3 )| 

y/3'lcr 


1 2 


(19) 


where SGi and SG2 are the outputs of the axial 
strain gages placed respectively on the lower and up- 
per outermost fibers of the tube at the panel center 
region, and RSG2 and RSG3 are the outputs of delta- 
rosette strain gage legs other than the leg parallel 
to the tubes. The bar indicates the initial nonzero 
strain gage readings at elevated temperatures when 
no mechanical loads were applied (panels were soaked 
at uniform temperature to determine these initial 
strains due to gage drift, apparent strain, and 
unintentional temperature nonuniformity). 

Panel Buckling Loads 

After the buckling load F cr is determined using 
the F/S method, the associated panel axial compres- 
sion stress resultant at buckling (N x ) C r and panel 
shear stress resultant at buckling {N xy ) cr must be 
determined. If N* and N* y are respectively the 
panel axial compression and shear stress resultants 
associated with the maximum applied load F* (see 
fig. 7 ), a nd if ( RSGf — RSGi), (RSG£ - RSG 2 ), and 
(RSG3 — RSG3) are the readings of the three legs of 
the rosette strain gage when F = F* ((RSGj— RSGi) 
being in the axial direction), then N* and N* y can 
be calculated as 


N* = Et{ RSG^-RSGi) (20) 


N* xy = ±Gt\(RSG * 2 - RSG 2 ) - (RSG3 - RSG 3 )| 

( 21 ) 

where t — 0.0368 in. is the equivalent extensional 
thickness of the panel, t = 0.016 in. is the thickness 
of the tubular wall, and RSG* (i = 1 to 3) are the 
rosette strain gage readings at F — F * . 

If the extrapolation factor k (see fig. 7) is defined 
as 



{N x ) cv = k N* 

(23) 

{N X y ) cr — k N* y 

(24) 


The values of ( N x ) cr and ( N xy ) cv thus obtained 
from test data are used in constructing the buckling 
interaction figures. 

Test Equipment 

Combined Loads 

Hypersonic wing test structure . The hypersonic 
wing test structure (HWTS), shown in figure 8, has 
a planform area of 85 ft 2 and is a portion of a 
proposed hypersonic research airplane (HRA) wing 
shown in figure 9. The HWTS was constructed 
based on the knowledge gained from the study of hot 
structural concepts for a Mach 8 hypersonic cruise 
vehicle with a 2.5g pull-up capability (refs. 1 and 
12). The HWTS was tested extensively in the past 
(ref. 13) to evaluate the hot-wing structural concept 
and to evaluate flight loads instrumentation, high- 
temperature calibration methods, and temperature 
simulation techniques. The beaded skin panels and 
corrugated spars and ribs are made of Rene 41, a 
nickel-base alloy. The heat shields are single-sheet 
panels which are slightly corrugated in the chordwise 
direction and are made of Rene 41 alloy except for 
those along the leading edge, which are made of 
TD Ni-20Cr. The Rene 41 heat shields are designed 
for locations where the surface temperature is less 
than 1800°F, and those made of TD Ni-20Cr are 
capable of operating with surface temperatures in 
excess of 1800°F. The heat shields are separated from 
the beaded skin panel by z-shaped support clips in 
order to minimize heat conduction from the heat 
shields to the substructure. The HWTS is connected 
to the support structure through a transition section 
and is mounted inverted so that wing loads produce 
compression on the lower surface of the HWTS. The 
transition section provides a load distribution buffer 
between the support structure and the test portion of 
the wing. The upper wing root zone (lower surface of 
HWTS) is the most highly compression-loaded area, 
and the five beaded panels there were replaced with 
five tubular panels for the nondestructive buckling 
tests. Figure 10 shows the HWTS with the heat 
shields removed to reveal the substructure and the 
z-shaped clips for supporting the heat shields. 

Mechanical loading system. Figure 11 shows the 
location of the applied load points on the HWTS and 
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the locations of the five test tubular panels. Twenty 
closed-loop channels were used to control electrohy- 
draulic equipment which applied mechanical loads to 
the test structure at the load points. Ten hydraulic 
jacks were used to apply vertical loads (simulation 
of lift load) to the HWTS to induce compression 
loads in the test panels. Eight of those jacks ap- 
plied loads through two-point whiffletrees. Horizon- 
tal loads (simulation of drag and thermal loads from 
adjacent vehicle structure) were applied with the re- 
maining 10 hydraulic jacks at single load points at 
the fore and aft edges of the HWTS to induce shear 
loads in the test panels. Pressure loads (which in- 
duced bending loads in the panels) were applied nor- 
mal to the upper surface of each test panel by using a 
0.003-in.-thick stainless-steel pressure pan positioned 
over each test panel. Each pan, which was bolted to 
the perimeter of a panel, thus formed one side of a 
pressure box. (See fig. 12.) 

Heating system . The system used to simulate 
aerodynamic heating of the HWTS is shown in fig- 
ure 13. The system was designed to heat the entire 
upper and lower surfaces of the HWTS to the tem- 
peratures corresponding to a Mach 8 flight profile. 
Infrared quartz lamps mounted on water-cooled pol- 
ished aluminum reflectors (as shown in fig. 14) were 
used to provide radiant heat. The system consisted 
of separate lower and upper heating units which were 
slightly contoured to match the surface shape of the 
HWTS. The units were mounted on rollers and tracks 
(see fig. 13) so that they could be easily removed 
for access to the HWTS and then be precisely repo- 
sitioned. The heating units were positioned with 
the reflector surfaces approximately 6 in. from the 
heat shields of the HWTS. Gaps in the lower heater 
were provided along the spar caps to allow clearance 
for connectors from the vertical loading system (see 
fig. 13). To fill in those gaps between load points, 
a double row of quartz lamps mounted on separate 
long, narrow water-cooled aluminum reflectors (i.e., 
strip heaters) was installed parallel to the spar caps. 
The temperatures of the panels were controlled by 
signals to the heating system from feedback thermo- 
couples attached to the heat shield exterior surfaces. 
The plumbing for the reflector cooling water included 
a pressure gage for each feed line to assure adequate 
coolant pressure. During the elevated-temperature 
tests, insulation curtains were draped around the 
HWTS and the heating system to reduce radiative 
and convective heat losses. (See fig. 13.) 

Compression 

A universal tension-compression testing machine 
was used for individual-panel axial compression 


buckling tests to obtain additional room-temperature 
buckling data in pure compression. Figure 15 shows 
the test machine with a test panel mounted. A to- 
tal of 11 displacement transducers (DT’s) were used 
to measure the out-of-plane deformations of the test 
panel. 

The surfaces of the upper and lower platens 
(which come into direct contact with the panel ends) 
were machined flat to ensure pure compression load- 
ing and to eliminate possible bending because of mis- 
aligned platen surfaces. The lower platen rested on 
a spherical seat and provided proper alignment with 
the test panel. 

End supports mounted on the panel provided 
surfaces for load transfer and served as reinforcement 
for the elimination of warping of the panel ends. The 
surfaces of both end supports were milled parallel 
with each other and perpendicular to the panel tube 
axes to provide pure compression load transfer. 

The panel vertical edges were bolted to the 
z-section stiffeners to approximate the stiffness con- 
ditions of a wing- mounted support. The interface 
between the panel and the stiffeners was lubricated. 
The holes on the stiffeners were oval shaped so that 
the bolts could move when the panel deformed. 

Strain Gage Instrumentation 

The strain gage locations on the surfaces of the 
five test panels are shown in figure 16. The view 
in the figure is looking downward from the top of 
the test panels. The strain gages with parentheses 
were located on the upper surfaces of the panels, and 
the rosette strain gages with square brackets were 
used for the elevated-temperature tests. The strain 
gages on the upper and lower surfaces of the panel 
tubes were single axial strain gages of two types: 
(1) foil type (circular symbol) and (2) capacitance 
type (square symbol). Of the axial strain gages, 
only the capacitance strain gages were capable of 
operating at temperatures above 550° F. The strain 
gages on the surface of the panel flat areas were the 
delta-rosette foil type and are indicated by the tri- 
angles in figure 16. At temperatures above 550° F, 
the bonded rosette gages were replaced with welded 
gages which are capable of operating at a tempera- 
ture of 1200°F. The delta-rosette strain gages were 
used to make measurements at three angular orienta- 
tions spaced 120° apart starting in the direction par- 
allel to the wing spars and rotating clockwise (when 
looking down on the test panels and inboard). The 
accuracy of the data acquisition system for strain 
gage measurements was ±5 x 10 -6 , which represents 
0.3 percent of the strain gage calibration output. Fig- 
ure 17 shows the full instrumentation of strain gages 
and thermocouples on test panel 5, and figure 18 
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shows the fully instrumented test panels attached to 
the HWTS lower wing root test area with panel 3 
removed to show the pressure pan interior. 

For the elevated-temperature tests, the strain 
gage outputs were corrected by subtracting the ini- 
tial nonzero readings at temperature without me- 
chanical load. Figure 19 shows the strain produced 
when the weldable gages are welded to Rene 41 and 
heated. This apparent strain would totally account 
for the initial nonzero reading if no gage drift or 
strain due to thermal stress exists. These initial read- 
ings were generally of the magnitude shown in fig- 
ure 19, indicating that gage drift and thermal stress 
were small. Figure 20 shows the full instrumentation 
of strain gages on the outer surface of test panel 1 for 
room-temperature, pure-compression, single-panel 
buckling tests. 

Test Procedure 

Combined Load Tests 

To generate a wide range of buckling data, a 
series of nondestructive buckling tests using the F /S 
method was conducted under various combined load 
conditions and at three temperatures (70°F (room 
temperature), 550° F, and 1000° F). Table 2 shows 
the maximum loads applied at the load points for 
different load conditions. Before the series of tests at 
70°F, the pressure system was checked to assure that 
a constant pressure level could be maintained during 
the tests. The pressure load was always maintained 
at the constant level of 0.75 psi or at 0 psi. Before the 
elevated-temperature tests, the heating system was 
checked to assure that constant temperature levels 
could be maintained over one test period. During 
the combined load, elevated-temperature tests, heat 
was first applied to raise the HWTS wing panels 
to a uniform temperature, and then pressure and 
mechanical loads were applied in that order. Table 3 
shows the test numbers and the corresponding load 
conditions. 


Data Reduction 


In applying the F/S method mentioned previ- 
ously, a typical vertical load (lift force) F703 located 
at load point 703 (associated with jack number 3; see 
fig. 11) was selected as F in equation (16) (or fig. 7) 
in the F/S calculations for all the test panels. The 
load F703 was arbitrarily selected as a representa- 
tive measurement of all applied loads since all loads 
were directly proportional to each other and were 
applied simultaneously. For F — F703, equation (22) 
becomes 


k = 



(^703)cr 

Z7* 

■*703 


(25) 


where (FVo3)cr is the value of F703 at the predicted 
buckling point and Fj 03 is the maximum value of 
F703 in the nondestructive buckling test. The buck- 
ling values of the panel stress resultants (N x ) cr and 
{N xy ) cr may be calculated by using equations (20), 
(21), (23), (24), and (25). For example, for panel 1 
at room temperature with strain gage combination 
RSG933, RSG934, and RSG935 (see fig. 16), equa- 
tions (20) and (21) may be written as 


iV* - Et[{ RSG* 933) - 0] (26) 


N* = — ^=G£|[(RSG* 934) - 0] - [(RSG* 935) - 0]| 


V3 


(27) 


and (N x ) cv and (N X y)er can be calculated from 
equations (23) and (24) as follows: 


(N x )cr = kN* = (28) 

■*703 

(N X y)cr = kN* y = (29) 

/03 


Results and Discussion 

Combined Loads 


Compression Tests 

Because the nondestructive buckling tests failed 
to produce results in pure compression at room tem- 
perature with p — 0 psi (see Results and Discussion ), 
two panels (panels 1 and 3) were tested to buckling 
failure at that load condition in a universal tension- 
compression testing machine. During the tests, the 
signals from the load cell, strain gage, and deflec- 
tometer channels were fed into the data acquisition 
system so that F/S plots could be generated. The 
buckling loads obtained from the F/S plots could 
then be compared with the actual buckling loads. 


Figures 21, 22, and 23 show respectively the 
force/stiffness (F/S) plots for the three typical 
tests 4.2.6 ( p — 0.75 psi), 4.4.6 (p = 0.75 psi), and 
4.3.4 (p = 0 psi) at 70°F. The solid curves shown 
in the figures were drawn from least-squares fits of 
the test data points. For some tests, or for certain 
strain gage combinations in the same test, the least- 
squares-fit curves based on equation (17) started to 
bend upward immediately after the last data points 
and intersected with the limit strain lines at points 
predicting rather high values of the buckling loads. 
For such occurrences, the extrapolations of the test 
data curves were accomplished by visually fitting the 
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data. No attempt was made to improve the extrap- 
olations by modifying equation (17) or by excluding 
data at low loads from the data set used to define 
the least-squares fit. The dashed curves shown in 
the three plots in figures 21(a), 21(c), and 21(d) were 
drawn as visual extrapolations. Notice that the plot 
in figure 21(b) shows excellent extrapolation of the 
least-squares-fit curve. With the existence of the lat- 
eral pressure (p = 0.75 psi), the F/S plots are usually 
convex upward. (See figs. 21 and 22.) However, when 
there is no lateral pressure (p = 0 psi), the F/S plots 
usually give strong convex downward curves except 
for the low-load region, giving quite accurate buck- 
ling load predictions (see fig. 23). The F/S plots 
for the rest of the tests where buckling loads are 
predicted are similar and therefore are not shown. 

Compression 

Figure 24 shows panel 3 after the room- 
temperature single-panel compression buckling test. 
A loud popping sound and a noticeable drop in load 
occurred at buckling. On one side of the panel, local 
failures (bead cripplings) were observed at three of 
the four beads within 5 in. of the panel center. (See 
fig. 24.) Only two small creases on a bead were ob- 
served on the other side. None of the strain gages 
were located directly at a buckle, although one small 
buckle not visible in figure 24 was located near strain 
gage 515. The behavior of panel 1 was similar except 
that buckles occurred in all four beads. 

Figure 25 shows out-of-plane displacements of 
three cross sections of test panel 3 at the panel com- 
pression load of ( N x )cr = 2138 lbf/in. immediately 
before buckling. The smaller displacements mea- 
sured near the panel edges indicate the existence 
of a stabilizing effect from the edge supports and 
may explain why the panel did not fail in general 
buckling at the predicted room-temperature Euler 
wide-column load of 1684 lbf/in. obtained from ref- 
erence 6. The local buckling load of 1622 lbf/in., 
also determined from reference 6, was exceeded by 
an even greater amount. As discussed subsequently, 
these results indicate that the theory used to design 
and analyze these panels appears to be unnecessarily 
conservative in compression. 

Figures 26 and 27 show respectively the F/S 
plots for panels 3 and 1 for the room-temperature 
single-panel compression tests. Notice that the ac- 
tual buckling points are located in the vicinity of 
the limit strain lines, and a visual extrapolation of 
the F/S data shown in figure 26(b), which were from 
strain gages located near a buckle, would give excel- 
lent agreement with the failure force of 41051 lbf. 
These results indicate that the F /S method could 


fairly accurately predict buckling loads in pure com- 
pression. (As discussed later, F/S predictions of 
buckling failure in pure compression at room tem- 
perature were not obtained for panels tested in the 
HWTS (fig. 28(a)) because the applied load was 
limited to a low value.) 

Comparison With Theory 

Table 4 summarizes the results of all the tests. In 
the table, N* and N* y are associated with F$ { q 3 and 
(N x ) cr and {N xy ) cr with (FVo3)cr- For most tests, 
the extrapolation factor k = (^703)^/^703 was ^ e " 
tween 2 and 3, thus indicating relatively large ex- 
trapolations. The large extrapolations were neces- 
sary because the applied loads were limited to less 
than 50 percent of the wing panel buckling load to 
prevent failure of the spar flanges, which were, by 
design, the critical components of the HWTS. 

The results given in table 4 were used to construct 
the buckling interaction plots shown in figure 28 
for different temperatures with or without pressure. 
The theoretical buckling interaction curves shown in 
the figures for comparison were generated through 
temperature and material modulus corrections of the 
results given in table 1 of reference 6. The two 
curves shown for local buckling were plotted with 
and without the additional knockdown factor of 0.86 
for equation (5). 

For high compression (see fig. 28(b)), the max- 
imum applied loads were not large enough to give 
accurate buckling data through the F/S data ex- 
trapolations. Nevertheless, the correlations between 
the test data and the predictions are fairly good 
in spite of data scatter resulting from the large ex- 
trapolations for some tests (e.g., near the iV x -axis). 
Most of the data points, including the actual buck- 
ling points obtained from the single-panel compres- 
sion tests (see fig. 28(a)), fall outside the predicted 
interaction curves, indicating the theoretical results 
are conservative. (As previously mentioned, these re- 
sults are based on a thickness of 0.016 in. and include 
effects from an assumed 0.001L initial imperfection). 
In all other tests, the theoretical curves fall within 
the scatter of the experimental data. Therefore, the 
experimental buckling data verify the theory for the 
applied test conditions and indicate that the addi- 
tional knockdown factor of 0.86 for local buckling 
in compression recommended in reference 6 is not 
necessary. 

As would be expected, the existence of lateral 
pressure, which adds a bending stress to the com- 
pression stress, decreases the compression buckling 
load ( N x ) cv considerably, but only slightly decreases 
the shear buckling load (N xy ) C r- Also, the buck- 
ling interaction curve shrinks as the temperature 
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is increased because of the decreases in E and G. 
Finally, the room-temperature pure-compression 
buckling loads per unit panel weight ( N x ) cr /w 
(w = 0.0161 lbm/in 2 ) for the tubular panels 1 and 
3 are 1.5067 x 10 5 in. and 1.3280 x 10 5 in., re- 
spectively. These values are slightly higher than 
(N x ) cr/w = 1 1507 X 10 5 in. (w = 0.0146 lbm/in 2 ) 
reported in reference 9 for the Rene 41 beaded panels 
which were originally used on the HWTS. 

Conclusions 

Five Rene 41 tubular panels which show promise 
of low structural mass and high buckling strength 
were installed as replacement root-chord wing panels 
on a section of a hot hypersonic wing test structure. 
To characterize their buckling behavior, the panels 
were exposed to nondestructive buckling tests under 
different combined load conditions and different tem- 
perature environments representative of those which 
would be encountered in a hot hypersonic wing, ex- 
cept that the structure was maintained at uniform 
temperatures. Thus, the. results included the effect 
of changes in modulus with temperature while the 
complex thermal stresses which can arise when tem- 
peratures are not uniform were minimized. Buckling 
loads for the wide range of loads and temperatures 
were obtained without failing the test panels through 
use of the force/stiffness method. 

In spite of some data scattering because of large 
extrapolations, the overall test data correlated fairly 
well with theoretically predicted buckling interaction 
curves. The existence of lateral pressure added a 
bending stress to the compression stress and thereby 
decreased the compression buckling load (N x ) cr con- 
siderably. However, the effect of the lateral pressure 
on the reduction of the shear buckling load ( N xy ) CY 
was quite small. Also, increasing the temperature 
decreased both (N x ) cv and (A^ xy ) cr because of re- 
ductions in the shear modulus and the modulus of 
elasticity at elevated temperatures. The fact that al- 
most all the test data for nearly pure compression 
at room temperature fell outside the predicted buck- 
ling interaction curves indicates that the theory used 
to design and predict the buckling of the panels is 
conservative for that condition. For all other test 
conditions, the force/stiffness test data verified the 
theory and showed that the structural efficiency of 
the tubular panel is slightly higher than that of the 
beaded panel which it replaced. 
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Table 1. Material Properties and Theoretical Local Buckling 
Strains for Rene 41 


Parameter 

Values for temperatures of— 


70° F 

550° F 

1000°F 

E, psi 

31.6 X 10 6 

29.1 x 10 (> 

26.5 x 10 e 

G, psi 

12.1 x 10 6 

11.2 x 10 6 

10.2 x 10 6 

fey, PSI 

125.0 X 10 3 

120.0 x 10 3 

117.0 x 10 3 

( € c)cr 

2.330 x 10- 3 (t73 = 1) 

2.330 x 10~ 3 (??3 = 1) 

2.330 x 10~ 3 (t73 = 1) 

( e 6)cr 

3.412 x 10- 3 (r? 3 = 0.96) 

3.412 x 10- 3 (r? 3 = 0.97) 

3.412 x Hr 3 (r /3 = 0.99) 

7cr 

2.208 x 10~ 3 (?72 = 1) 

2.208 x 10~ 3 (??2 = 1) 

2.208 x 10 _3 (?72 = 1) 


Table 2. Jack Loads Applied to HWTS for Different Load Conditions 

(a) T — 70°F 


Jack 

Jack position 


Maximum load, a 

lbf, for load condition — 


3.1 

3.6 

3.8 

4.1 

4.2 

4.3 

4.4 

i 

Vertical 

1680 

2521 

4 201 

3500 

2800 

1200 

-4000 

2 

Vertical 

867 

1304 

2173 

6000 

2800 

2000 

2400 

3 

Vertical 

624 

735 

1557 

6000 

2800 

2000 

2400 

4 

Vertical 

863 

1274 

2157 

3500 

2800 

-2800 

-4000 

5 

Horizontal 


1650 

2 750 


-5000 

-5000 

-5000 

6 

Horizontal 


4260 

7100 


-6500 

-6500 

-6500 

7 

Horizontal 


-4030 

-6 716 


6500 

6500 

6500 - 

8 

Horizontal 


-1232 

-2054 


6500 

6500 

6500 

9 

Horizontal 


-3840 

-6350 


6000 

6000 

6000 

10 

Horizontal 


731 

1227 


6000 

6000 

6000 

11 

Vertical 

2265 

3398 

5663 

3000 

800 

-800 

-4000 

12 

Vertical 

-95 

-193 

-238 

1800 

1400 

1200 

1200 

13 

Horizontal 


-6692 

-11 154 


6500 

6500 

6500 

14 

Horizontal 


-3945 

-6575 


6500 

6500 

6500 

15 

Horizontal 


-1740 

-2900 


-6500 

-6500 

-6500 

16 

Horizontal 


2400 

4000 


-6500 

-6500 

-6500 

17 

Vertical 

1267 

1901 

3168 

3500 

2800 

2800 

4000 

18 

Vertical 

-1896 

-2770 

-4616 

3500 

4800 

4800 

4000 

19 

Vertical 

1760 

3640 

4400 

2800 

2800 

2500 

4000 

20 

Vertical 

592 

814 

1356 

1800 

1400 

1200 

1200 


a Positive values indicate tension; negative values indicate compression. 
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Table 2. Continued 
(b) T = 550°F 


Jack 

Jack position 

Maximum load, a lbf, 

:or load condition— 

6.6 

6.8 

7.1 

7.2 

7.3 

7.4 

i 

Vertical 

2 521 

4201 

3500 

2800 

1200 

-4000 

2 

Vertical 

1304 

2173 

6000 

2800 

2000 

2900 

3 

Vertical 

735 

1559 

6000 

2800 

2000 

2400 

4 

Vertical 

1294 

2157 

3500 

2800 

-2800 

-4000 

5 

Horizontal 

1650 

2 750 


-5000 

-5000 

-5000 

6 

Horizontal 

4260 

7100 


-6500 

-6500 

-6500 

7 

Horizontal 

-11030 

-6716 


6500 

6500 

6500 

8 

Horizontal 

-1232 

-2054 


6500 

6500 

6500 

9 

Horizontal 

-3840 

-6350 


6000 

6000 

6000 

10 

Horizontal 

737 

1227 


6000 

6000 

6000 

11 

Vertical 

3398 

5663 

3000 

800 

-800 

-4000 

12 

Vertical 

-143 

-238 

1800 

1400 

1200 

1200 

13 

Horizontal 

-6692 

-11154 


6500 

6500 

6500 

14 

Horizontal 

-3945 

-6575 


6500 

6500 

6500 

15 

Horizontal 

-1740 

-2 700 


-6500 

-6500 

-6500 

16 

Horizontal 

2 400 

4000 


-6500 

-6500 

-6500 

17 

Vertical 

1901 

3168 

3300 

2800 

2800 

4000 

18 

Vertical 

-2 770 

-4616 

3500 

4800 

4800 

4000 

19 

1 Vertical 

2640 

4 400 

2800 

2800 

2800 

4000 

20 

Vertical 

814 

1356 

1800 

1400 

1200 

1200 


“Positive values indicate tension; negative values indicate compression. 




Table 2. Concluded 
(c) T - 1000° F 


Jack 

Jack position 


Maximum load,® lbf, 

for load condition — 


9.3 

9.5 

8.1 

8.2 

8.3 

8.4 

1 

Vertical 

2 521 

4 201 

3500 


1200 

mssm 

2 

Vertical 

1304 

2173 

6000 


2000 

■SI 

3 

Vertical 

935 

1559 

6000 


2000 

B9 

4 

Vertical 

1294 

2157 

3500 

2800 

-2800 

-4000 

5 

Horizontal 

1650 

2 750 


-5000 

-5000 

-5000 

6 

Horizontal 

4 260 

7100 


-6500 

-6500 

-6500 

7 

Horizontal 

-4 030 

-6 716 


6500 

6500 

6500 

8 

Horizontal 

-1232 

-2054 


6500 

6500 

6500 

9 

Horizontal 

-3 840 

-6 350 


6000 

6000 

6000 

10 

Horizontal 

737 

1229 


6000 

6000 

6000 

11 

Vertical 

3398 

5 663 

3000 

800 

-800 

-4000 

12 

Vertical 

-143 

-238 

1800 

1400 

1200 

1200 

13 

Horizontal 

-6672 

-11154 


6500 

6500 

6500 

14 

Horizontal 

-3945 

-6575 


6500 

6500 

6500 

15 

Horizontal 

-1740 

-2 700 


-6500 

-6500 

-6500 

16 

Horizontal 

2 400 

4 000 


-6500 

-6500 

-6500 

17 

Vertical 

1901 

3168 

3500 

2800 

2800 

4000 

18 

Vertical 

-2 770 

-4616 

3500 

4800 

4800 

4000 

19 

Vertical 

2 640 

4 400 

2800 

2800 

2800 

4000 

20 

Vertical 

814 

1356 

1800 

1400 

1200 

1200 


“Positive values indicate tension; negative values indicate compression. 



Table 3. Test Numbers and 


Test 

Load 

condition 

4.1.4 

4.1 

4.2.6 

4.2 

4.3.3 

4.3 

4.4.6 

4.4 

4.2.7 

4.2 

4.3.4 

4.3 

4.4.7 

4.4 

6.8.3 

3.8 

7.2.4 

4.2 

7.3.5 

4.3 

7.4.4 

4.4 

7.1.8 

4.1 

8.2.2 

4.2 

8.2.2 

4.2 

8.3.2 

4.3 

8.3.2 

4.3 

8.3.5 

4.3 

8.3.5 

4.3 

8.4.6 

4.4 

8.4.6 

4.4 

8.1.3 

4.1 

8.1.3 

4.1 

Single-panel test 1.1 

(a) 

Single-panel test 3.1 

(a) 


a N x 7 ^ 0 ; N X y — 0 . 


Load Conditions 


Pressure, 

psi 

Temperature, 

°F 

0.75 

70 

.75 

70 

.75 

70 

.75 

70 

0 

70 

0 

70 

0 

70 

.75 

550 

.75 

550 

.75 

550 

.75 

550 

.75 

550 

0 

1000 

.75 

1000 

0 

1000 

.75 

1000 

0 

1000 

.75 

1000 

0 

1000 

.75 

1000 

0 

1000 

.75 

1000 

0 

70 

0 

70 
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Table 4. Summary of Test Results 
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a Location 1 - half-panel location; location 2 — quarter-panel location. 
b From visual extrapolation of F/S data. 
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Figure 2. Rene 41 tubular panel with heat-shield supports removed. 
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Figure 3. Plasticity correction curves for local buckling of Rene 41 circular arc element for T = 70° F. 
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Figure 4. Plasticity correction curves for local buckling of Rene 41 circular arc element for T = 550°F. 
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Figure 5. Plasticity correction curves for local buckling of Rene 41 circular arc element for T — 1000°F. 
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Region of At buckling failure 



Figure 7. Force/stiffness plot for local buckling. 
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igure 8. Hypersonic wing test structure (HWTS). 
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Figure 9. Hypersonic wing test structure as part of hypersonic research airplane. 
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gure 10. Hypersonic wing test structure with heat shields removed to show substructure. 
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Figure 11. Applied load distribution on HWTS and locations of five test tubular panels. 
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Figure 13. Hypersonic wing test structure combined mechanical and thermal loading test setup. 
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Figure 15. Tubular panel installed in testing machine for axial compression buckling test. 
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Figure 18. Four Rene 41 tubular panels attached to hypersonic wing test structure for buckling tests. View 
looking up and aft at the lower side of test structure. Panel 3 removed to show pressure pan interior. 
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Figure 21. Force/stifFness plots for panel 3 test 4.2.6. T = 70°F; p = 0.75 psi. 
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Figure 22. Force/stiffness plots for panel 3 test 4.4.6. T = 70°F; p = 0.75 psi. 
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1% = 1-22- (b) = 1.22. 



Figure 23. Force/stifFness plots for panel 3 test 4.3.4. T = 70°F; p = 0 psi. 
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Figure 24. Buckled panel 3 after compression buckling test. 
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Figure 25. Out-of-plane displacements of test panel 3 immediately before buckling for single-panel compression 
test. 
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(b) T = 70°F; p = 0.75 psi. Strain gage set refers to locations 1 and 2 of table 4. 

Figure 28. Continued. 
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Figure 28. Continued. 
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